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ABSTRACT 

Recent high-resolution measurements suggest that the soft X-ray spectrum of 
obscured Radio Galaxies (RG) exhibits signatures of photoionised gas (e.g. 3C 445 
and 3C 33) similar to those observed in radio-quiet obscured Active Galactic Nuclei 
( AGN) . While signatures of warm absorbing gas covering a wide range of temperature 
and ionisation states have been detected in about one half of the population of nearby 
Seyfert 1 galaxies, no traces of warm absorber gas have been reported to date in the 
high-resolution spectra of Broad Line Radio Galaxies (BLRG). We present here the 
first detection of a soft X-ray warm absorber in the powerful FRII BLRG 3C 382 using 
the Reflection Grating Spectrometer (RGS) on-board XMM-Newton. The absorption 
gas appears to be highly ionised, with column density of the order of 10 22 cm -2 , 
ionisation parameter log£ >2 erg cm s _1 and outflow velocities of the order of 10 3 km 
s _1 . The absorption lines may come from regions located outside the torus, however at 
distances less than 60 pc. This result may indicate that a plasma ejected at velocities 
near the speed of light and a photoionised gas with slower outflow velocities can coexist 
in the same source beyond the Broad Line Regions. 
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1 INTRODUCTION 

The launch of XMM-Newton and Chandra both carrying on- 
board diffraction grating spectrometers with high resolving 
power, dramatically improved our knowledge of the gaseous 
environment of Active Galactic Nuclei. At first, very bright 
and nearby Seyfert galaxies were explored and many results 
obtained in the last decade. Depending on the observer's 
line-of-sight, different features could be detected. In X-ray 
obscured Seyfert 2s high-resolution observations have re- 
vealed soft X-ray spectra dominated by line emission from 
H-like and He-like elements as Neon, Oxygen, Carbon, as 
well as from Fe L-shell transitions (Fe XVII to Fe XXIV) 
(Sako et al. 2000; Kinkhabwala et al. 2002; Bianchi et al. 
2005; Pounds & Vaughan 2006; Guainazzi & Bianchi 2007). 
The soft X-ray emission is most likely produced by photoion- 
isation in a gas located in the Narrow Line Region (NLR) 
(Bianchi, Guainazzi & Chiaberge 2006). 

Signatures of outflowing ionised gas have been detected 
in the X-ray spectra of about 50% of nearby Seyfert Is. 
The data indicate the presence of a multi-phase gas with 



an ionisation parameter £ * spanning three orders of magni- 
tude and column densities ranging between io 20-22 cm~ 2 , 
with velocities of the order of a few hundreds to thousands 
km s^ 1 (Blustin et al. 2005; McKernan et al. 2007). Where 
these winds originate and are powered is still an open ques- 
tion. Several regions at different distances from the central 
engine, like the NLR, the obscuring torus and the accretion 
disk have been proposed as a possible location for the gas 
(Crenshaw et al. 2003 and references therein). 

For Radio Galaxies the investigation of the nuclear 
environment through high-resolution spectroscopy has just 
started. There are some hints that in obscured RGs (seen at 
high inclination angles) the soft X-ray spectra are similar to 
their Radio-Quiet (RQ) counterparts, Seyfert 2s. Evidence 
for photoionised gas has already been revealed by emission 

, L 

£= -, L is the 1-1000 Rydbcrg source luminosity, n e is the 

n e R 2 

gas density and R is the distance from the source (Tarter Tucker 
and Salpeter 1969) 
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lines in the soft X-ray spectra of 3C 445 (Grandi et al. 2007; 
Sambruna et al. 2007) and 3C 33 (Torresi et al. 2009). In 
the latter case, the morphological coincidence between soft 
X-ray and optical [OIII]A5007 emission suggests the NLR as 
the location of such gas (Torresi et al. 2009). 

The current picture is less clear for the radio-loud (RL) 
counterparts of Seyfert Is, where Doppler amplification of 
the jet flux can play a key role in blurring the spectral fea- 
tures: for example, RGs are known to have weaker Compton 
reflection signatures and weaker Fe Kq lines when compared 
to Seyfert Is (Grandi, Malaguti & Fiocchi 2006 and refer- 
ences therein) . Up to now, the search for warm absorbers in 
unobscured RGs (seen at low inclination angles) has been 
unsuccessful, except for a claim of its existence in 3C 382 
and 3C 390.3 from analysis of ASCA data (Reynolds 1997). 
Therefore it is completely unclear whether winds can form 
and/or survive in the nuclear regions of RL AGNs, where 
large amounts of the gravitational power is funneled into 
relativistic collimated plasma. Addressing this point is fun- 
damental to clarify the physical link among accretion flows, 
winds and jet, and the influence that each of these compo- 
nents can have on the host galaxy. 

3C 382 is a nearby Radio Galaxy (z=0.0579), with lobe- 
dominated Fanaroff-Riley II radio morphology. It exhibits a 
long jet extending north-east of the core for 1.68' and two 
radio lobes, with a total extension of 3' (Black et al. 1992). 
At optical wavelengths 3C 382 shows strong and broad lines 
(FWZI>25,000 km s" 1 ), which are variable on timescales of 
months to years. An Hubble Space Telescope (HST) WFPC2 
image showed that this source is an elliptical galaxy strongly 
dominated by an unresolved nucleus (Martel et al. 1999). 
The optical emission of this source is clearly photoionised 
as can be inferred from the de-reddened line ratios, i.e. 
log[OIII]/H/?=0.90 and log[NII]/Ha=0.17 (Buttiglione et al. 
2009), compared with the diagnostic diagrams of Miller et 
al. (2003). 

In the X-ray band 3C 382 is a bright source (F2-iofc e v ~ 
3xl0 -11 erg cm -2 s _1 , from this XMM-Newton observa- 
tion). The 2-10 keV spectrum is well fitted with a single 
power-law; when the fit is extrapolated to lower energies 
a strong soft-excess is observed (Prieto 2000; Grandi et al. 
2001). ROSAT/HRI observations revealed extended X-ray 
(0.2-2.4 keV) emission around 3C 382 (Prieto 2000), that 
was confirmed by the spatial resolved study of Gliozzi et al. 
(2007). However, as discussed by Grandi et al. (2001) and 
Gliozzi et al. (2007) the extended thermal emission cannot 
account for the observed soft-excess, which must be instead 
of nuclear origin. 

Here we report the first high-resolution detection of an X-ray 
warm absorber in this BLRG. 



2 OBSERVATIONS AND DATA REDUCTION 

3C 382 was observed by XMM-Newton on 8th April 2008 
(OBSID: 0506120101) for about 40 ks. The data from EPIC 
(Striider et al. 2001), RGS (den Herder et al. 2001) and Op- 
tical Monitor (OM, Mason et al. 2001) instruments were op- 
crating. As we are primarily interested in the warm absorber 
features, only the RGS data (6-38 A) are presented in this 
paper. Data were processed with the task rgsproc of the SAS 
8.0.0, which combines the event lists from all RGS CCD, 



produces source and background spectra using standard ex- 
traction regions, and generates response matrices. The back- 
ground is estimated taking into account events from a region 
spatially offset from the source. All spectral fits presented 
in this work include absorption due to a line-of-sight Galac- 
tic column density of N//=7.0xl0 20 cm" 2 (Kalberla et al. 
2005). The cosmological parameters used throughout the 
paper are H = 70 km s _1 Mpc -1 , Q m = 0.3, = 0.7 
(Spergel et al. 2007). The spectral analysis was performed 
using the XSPEC v. 11. 3. 2 package (Arnaud 1996). 



3 RGS DATA ANALYSIS 

The first-order unbinned spectra from the two RGS cameras 
were simultaneously fitted in XSPEC using the C-statistic 
(Cash 1979). The quoted errors correspond to 90% confi- 
dence level (AC=2.71) for one interesting parameter. 

The 6-35 A (0.35-2 keV) continuum is well fitted by a 
single absorbed power-law with Y = 2.25±0.18 (C=5264 for 
4636 d.o.f.). A quite steep power-law is not unexpected in 
this source known to have a strong soft excess. We investi- 
gated whether an extended thermal component (APEC) is 
required by the data, as suggested by the Chandra analysis. 
The fit did not improve significantly, nor did the shape of 
the residuals (C=5261). Therefore we decided not to include 
any other continuum components. However, as discussed be- 
low, we checked for the presence of weak collisional emission 
lines. 

To start with a phenomenological approach, the spec- 
tra were divided into a number of regions containing about 
the same number of bins each (100). We performed an ac- 
curate inspection of every single region. Where the spectra 
showed a particularly evident absorption feature, a Gaus- 
sian profile with 0-width and a negative normalization was 
added to the model. Only lines yielding an improvement in 
the fit higher than AC>5 (corresponding to > 95% for two 
interesting parameters) were taken into account. For these 
structures the equivalent width was estimated and the corre- 
sponding atomic transition identified through several atomic 
databases such as CHIANTI (Dere et al. 1997; Landi et al. 
2006), ATOMDB WebGUIDE t and NIST *. The parame- 
ters of the most prominent lines are reported in Table 1 and 
the phenomenological fit is shown in Figure 1. Two other ab- 
sorption lines probably corresponding to FeXXI/MgXI and 
SIV/NVII transitions, although observed, are not reported 
in Table 1 because of their marginal detection (AC =5). 
None of the detected features is due to local (Galactic) ab- 
sorption. 

Finally we clearly detected an emission line (AC=6) 
identified as OVIII Lya (see Table 1). This is not an un- 
usual result. Emission features are often observed in sources 
showing signatures of warm absorber (Blustin et al. 2002; 
Turner et al. 2004). 

Once the presence of absorption lines was ascertained 
via the phenomenological analysis, the next step was to 
study the physical properties of the gas using a self consis- 
tent modeling of the spectrum. Extensive simulations were 

t http://cxc.harvard.edu/atomdb/WebGUIDE/index.html 
t http://physics.nist.gov 
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Table 1. Soft X-ray absorption and emission lines detected in the RGS spectrum of 3C 382. Columns: (1) Identification of the transition; 
(2) 3C 382 rest-frame wavelength (in A); (3) Laboratory wavelength (in A); (4) Emission-line intensity; (5) Equivalent width (in mA); 
(6) AC. 



Transition 


A(rest) 

(A) 


A (lab) 

(A) 


Flux 

(10 _5 photons cm -2 s _1 ) 


EW<°' 
(mA) 


AC M 


NeX 


12.07 


12.13 




°'-30 


AC=7 


FeXX 


12.78 


12.82 




4 o+3i 

*°-30 


AC=8 


OVIII Lya 


18.90 


18.97 




04+12 


AC=11 


OVIII Lya (em) 


18.99 


18.97 


^- ,u -3.32 


nq+16 
" zo -18 


AC=6 



(a) Errors are at the 90% confidence level for one interesting parameter. 

(b) Improvement in the C-statistic (Cash 1979) for adding a Gaussian line with two 
free parameters to a power-law model. 



carried out using the photoionisation code XSTAR (v2.11n9) 

§ . The output was a grid of models (afterwards simply called 
XSTAR model) that can be fitted to the data in XSPEC. 
Each model in the grid represents a photoionised spectrum 
characterized by a particular value of the column density 
and the ionisation parameter. Our explored parameter range 
is N H =[10 18 , 10 23 j and log£=[-3, 3]. We assumed a simple 
spherical geometry for the illuminated gas (covering factor 
= 1), solar abundances (Grevesse & Sauval 1998) , and an il- 
luminating continuum from 1 to 1000 Rydbergs derived from 
the intrinsic spectral energy distribution SED of the source, 
with L ion — 3 x 10 45 erg s _1 (Figure 2). The turbulent ve- 
locity of the gas was set to vt u , i,=100 km s _1 , as the limited 
RGS statistics did not allow to constrain the intrinsic width 
of the absorption features. 

The XSTAR model obtained from these simulations 
was applied as an XSPEC multiplicative component to the 
power-law and fitted to the data. The fit significantly im- 
proved (AC=160), as XSTAR can take into account any 
possible absorption feature in the spectrum. We note that 
the photon index is completely consistent with the one ob- 
tained from a power-law only fit. The presence of a highly 
ionised gas phase is confirmed, indeed the best fit required a 
high ionisation parameter log£;=2.69to'o3 er S cm s_1 an< i a 
column density N_f/=3.19to'79 xlO 22 cm -2 . The robustness 
of these results is attested by the N_ff- log£ confidence con- 
tours (68% 90% 99%) shown in Figure 3 (upper panel). We 
attempted to estimate the gas outflow velocities by measur- 
ing the line shifts between theoretical and observed values of 
the energy centroid allowing the redshift parameter (z) of the 
XSTAR model to vary. We found reasonable outflow speeds 
of the order of -1200l2oo km s _1 , consistent with typical 
outflow velocities in Seyfert Is (Blustin et al. 2002; Blustin 
et al. 2005) when the absorption gas is seen along the line- 
of-sight. Incidentally we point out that the phenomenogical 
identifications proposed in Table 1 are confirmed by the XS- 
TAR model (see Figure 3 lower panel). 

We note that the energy centroid of the OVIII Lya 



http: / /heasarc. nasa.gov/lheasoft/xstar/xstar. html 




Figure 1. Zoom of the 3C 382 RGS observed-frame spectra 
(black: RGS1; red: RGS2), in wavelength ranges, around the de- 
tected absorption and emission lines. The identified lines are: 
NeX-FeXX (upper panel), OVIII Lya in emission and absorption 
(lower panel). The green line is the best fit phenomenological 
model (Table 1). 



emission line does not require any blueshift (although with 
large uncertainties). However it may suggest that, contrary 
to the absorption lines observed along the line-of-sight, the 
emission comes from an extended region around the source, 
maybe a shell, expanding in all directions. 
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Figure 2. Radio to X-ray spectral energy distribution of 3C 382. 
Data are from literature (Rudnick et al. 1986; Knapp et al. 1990; 
Klavel et al. 2000; Lilly & Longair 1982; Chiaberge et al. 2000; 
Chiaberge et al. 2002; Grandi et al. 2001; Quillen et al. 2003; 
Ramos Almeida et al. 2007; Hinshaw et al. 2007 ). Optical and UV 
measurements de-reddened with the extinction curve of Cardelli 
et al 1989. The visual extinction is deduced by the Galactic col- 
umn density, assuming the gas-to-dust ratio from Shull et al. 1985 




UA) 

Figure 3. Upper panel: Contour plots for column density and 
ionisation parameter values of the warm absorber as obtained 
with the XSTAR model. The curves define the 68%, 90% and 
99% confidence intervals for two interesting parameters. Lower 
panel: Best fit rest-frame XSTAR model of 3C 382 in the overall 
RGS band (6-35 A) corresponding to N H =3.19l ( ! ) !,g xlO 22 cm" 2 
and log§=2.69±°'°3 erg cm s" 1 . 



4 DISCUSSION AND CONCLUSIONS 



The high-resolution analysis of the BLRG 3C 382 has re- 
vealed that relativistic collimated plasma and slower out- 
flows of photoionised gas can coexist in the same object. 
The outflowing gas is found in a hot phase as attested by 
the derived ionisation parameter log£=2.69 erg cm s -1 de- 
duced with the XSTAR simulations. Most of the light ele- 
ments are indeed completely ionised, and the strongest ab- 
sorption features correspond to H-like elements and L-shell 
iron states. In Seyfert Is warm absorbers are complex and 
multi-phase, often indicating that layers of gas at different 
ionisation states are part of the same outflow (Blustin et 
al. 2005). In the present case, only one component of hot 
gas was observed, without apparent traces of a colder com- 
ponent. Nonetheless, the physical conditions of the gas in 
3C 382 warm absorber are quite similar to those measured 
in the highest ionised absorber in their radio-quiet counter- 
parts (NGC 4051, Krongold et al. 2007; Mkr 279, Costantini 
et al. 2007; IRAS 13349+2438, Sako et al. 2001). Taking into 
account the short exposure time and the very simple as- 
sumptions in the XSTAR model, the result appears encour- 
aging. We are aware that setting v tU rb to an a priori fixed 
value could potentially overpredict the column density if the 
real turbulent velocity is larger. Interestingly, also the EWs 
of the absorption lines detected in 3C 382 are consistent with 
those generally measured in Seyfert Is (Yaqoob et al. 2003; 
Turner et al. 2004). This similarity suggests a negligible jet 
contribution to the soft X-ray spectrum, although the large 
uncertainties in EW prevent more quantitative statements 
from being made. 

In order to localize the warm absorber we estimate 
a minimum and a maximum distance for the gas follow- 
ing Blustin et al. (2005). The minimum distance is calcu- 
lated with the assumption that the outflow velocity has 

, , r, 2GM 

to be greater or equal than the escape one, R > — = — . 

"out 

For a M SH =l.lxl0 9 M Q (Marchesini et al. 2004), and av- 
erage outflow velocities of the order of 10 3 km s _1 (see 
Table 1) we deduce that the gas should be located at a 
distance larger than 10 pc. On the other hand, the max- 
imum distance of the gas can be obtained by combining 
the ionisation parameter and the column density. If the 
gas is confined in a thin layer (AR/R < 1 with AR the 
layer depth), and the volume filling factor equal to 1, it 

turns out that R < ( — ) < 60 pc. According to Ghis- 

ellini & Tavecchio (2008) the BLR and torus radii are ex- 
pected to scale with the square root of the disk luminosity: 
Rblr = W 17 L^ k 45 cm and Rtoms = 2.5 x 10 18 L^ 45 cm 
(Cleary et al. 2007). Thus for 3C 382 we can estimate typical 
sizes for the BLR and dusty torus region, i.e. Rblr —0.06 pc 
and Rtorus ~1.5 pc if Ldisk ~ Li 0n . These results suggest 
that the bulk of the highly ionised absorption gas, detected 
in 3C 382, is located beyond the torus, most likely in the 
NLR. Interestingly, we note that this is also the probable 
location of the photoionised emitting gas observed in 3C 33 
(Torresi et al. 2009), another powerful RG seen at higher 
inclination angles. 
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